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ABSTRACT 


The  inviscid,  linear,  non-axisymmetric,  temporal  stability  of  the  boundary  layer  associ¬ 
ated  with  the  supersonic  flow  past  axisymmetric  bodies  (with  particular  emphasis  on  long 
thin,  straight  circular  cylinders),  subject  to  heated  or  cooled  wall  conditions  is  investigated. 
The  eigenvalue  problem  is  computed  in  some  detail  for  a  particular  Mach  number  of  3.8, 
revealing  that  the  effect  of  curvature  and  the  choice  of  wall  conditions  both  have  a  significant 
effect  on  the  stability  of  the  flow. 

Both  the  asymptotic,  large  azimuthal  wavenumber  solution  and  the  asymptotic,  far  down¬ 
stream  solution  are  obtained  for  the  stability  analysis  and  compared  with  numerical  results. 
Additionally,  asymptotic  analyses  valid  for  large  radii  of  curvature  with  cooled/heated  wall 
conditions,  are  presented.  We  find,  in  general,  important  differences  exist  between  the  wall 
temperature  conditions  imposed  in  this  paper  and  the  adiabatic  wall  conditions  considered 
previously. 
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§1.  Introduction  and  Motivation 

The  recent  resurgence  of  interest  in  high  speed  flight  vehicles  has 
rekindled  interest  in  supersonic  and  hypersonic  flows.  One  important 
area  of  aerodynamic  study  is  that  of  boundary- layer  stability  and  the 
possible  transition  to  turbulence.  Since  turbulent  flows  result  in 
considerably  greater  skin  frictions  and  heat  transfer  coefficients  than  do 
laminar  flows,  any  method  by  which  a  boundary  layer  may  be  stabilized 
is  worthy  of  investigation. 

Generally  it  is  found  in  the  case  of  supersonic  boundary- layer  flows 
that  inviscid  disturbances  are  more  important  (i.e.  more  unstable)  than 
viscous  disturbances.  (Here,  we  characterise  inviscid  disturbances  as  being 
those  with  wavelengths  comparable  to  the  boundary- layer  thickness,  whilst 
viscous  disturbances  possess  much  longer  wavelengths).  This  is  in  contrast 
to  the  situation  encountered  for  many  incompressible  boundary- layer  flows 
where  viscous  instabilities  are  generally  dominant.  The  first  authors 
to  give  any  form  of  rigorous  mathematical  account  of  the  stability  of 
compressible  boundary  layers  were  Lees  and  Lin  (1).  Making  use  of 
a  rational  asymptotic  approximation  they  determined  that  the  quantity 

*  [  P  1  (where  u*  denotes  velocity  tangential  to  the  surface, 

dy  dy 

y*  the  normal  to  the  surface,  and  p*  the  fluid  density)  plays  a  key 
role,  similar  to  that  the  quantity  d^u*/dy*2  has  in  incompressible 
theory.  In  particular  at  the  point  where  the  above  expression  is  zero 
(y*=y i  ).  termed  the  generalised  inflexion  point,  then  there  may  exist  a 
neutral  mode  with  wavespeed  u  * ( y  j  * ) ;  neutral  modes  are  classed  as  being 
'subsonic  ' .  ''sonic''  or  'supersonic''  depending  on  how  the  freestream 
Mach  number  is  related  to  the  wavespeed  (Mack  (2)),  If  the  neutral 
disturbance  is  subsonic  then  the  mode  decays  in  the  far-fieid;  supersonic, 
neutral  disturbance  modes  exhibit  an  oscillatory  behaviour  in  the  far-field; 
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a  sonic  mode  occurs  at  the  crossover  point  between  subsonic  and  supersonic 
eases.  Mathematically,  these  classifications  are  directly  related  to  the 
non-dimensional  wavespecd  c  (defined  in  (2.10)  below),  and  the 
frcc-strcam  Mach  number  M^.  For  subsonic  disturbances  we  have 

1  '  *  c  <  1+  1 


for  sonic  disturbances  we  have 


c 


or 


c 


and  for  supersonic  disturbances  we  have 


c  < 


1- 


c  >  1  + 


1 

5L- 


Arguments  relating  to  generalised  inflexion  points  have  no  implications 
for  supersonic  neutral  modes. 

Lees  (3)  considered  the  effect  that  wall  cooling  has  on  the 
stability  of  compressible  boundary  layers  on  the  basis  of  asymptotic 
theory.  He  predicated  that  with  sufficient  wall  cooling  the  boundary 
layer  could  be  completely  stabilized  and  presented  a  criterion  whereby 
the  ratio  of  wall  temperature  to  the  recovery  temperature  at  which  the 
critical  Reynolds  number  becomes  infinite,  can  be  computed.  Even  though 
Lee’s  original  work  contained  numerical  errors,  subsequent  authors 
including  Van  Driest  (4)  and  Dunn  and  Lin  (5)  showed  that  Lee’s  predictions 
appeared  to  be  correct. 

Van  Driest  (4)  calculated  the  cooling  required  to  completely 
stabilize  the  boundary  layer  on  a  flat  plate  at  supersonic  speeds  with 
zero  pressure  gradient.  Whereas  Lee’s  investigations  were  limited  to  low 
supersonic  flows.  Van  Driest  predicated  that  complete  stabilization  was 
achievable  by  wall  cooling  over  a  wide  range  of  Mach  numbers  up  to 
hypersonic  flow.  He  found,  however,  that  for  Mach  numbers  greater  than  9, 
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it  is  impossible  to  stabilize  the  boundary  layer  with  any  amount  of 
cooling  when  a  Prandtl  number  of  0.75  and  the  Sutherland 
vi scosi ty- temperature  law  are  assumed. 

The  above  predictions  of  Lees  (3)  and  Van  Driest  (4)  were  based 
on  the  asymptotic  theory  of  two-dimensional  disturbances.  Dunn  and  Lin 
(5)  extended  this  work  to  include  three-dimensional  disturbances  and 
they  indicated  that  cooling  was  indeed  an  effective  method  by  which  the 
boundary  layer  could  be  stabilized  for  moderate  supersonic  Mach  numbers. 
Based  on  their  asymptotic  analysis,  Dunn  and  Lin,  however  (wrongly) 
concluded  that  at  supersonic  free-stream  Mach  numbers  the  boundary  layer 
can  never  be  completely  stabilized  by  moling  with  respect  to  all  three- 
dimensional  disturbances. 

However  the  above  asymptotic  work  is  found  to  have  its  limitations. 

The  theory  above  indicates  the  possible  existence  of  only  subsonic 
generalised  modes  for  cooled  wall  conditions.  In  a  series  of  papers  Mack 
(2,  6.  7,  8,  9,  10,  11),  using  numerical  techniques,  demonstrated  that 
for  compressible  boundary  layers  there  in  actual  fact  exist  a  large 
(probably  infinite)  number  of  unstable  modes,  including  many  with  supersonic 
neutral  points,  which  may  not  be  predicated  by  the  above  asymptotic  work. 

For  the  case  of  three-dimensional  disturbances  in  Blasius-type  boundary 
layers.  Mack  showed  that  the  so-called  ’’first  mode’’  of  instability 
undergoes  stabilization  with  sufficient  wall  cooling,  thus  verifying  the 
above  mentioned  asymptotic  analysis.  However,  Mack  found  in  the  case  of  the 
’’second  mode’’  of  instability  (which  turns  out  to  be  more  important  because 
of  the  larger  growth  rates  involved)  cooling  the  wall  actually  destabilizes 
this  mode.  According  to  recently  presented  numerical  observations  by  Mack 
(2,  10)  oblique  disturbances  can  be  completely  stabilized  by  wall  cooling 
for  Mach  numbers  up  to  5.8,  although  they  do  require  a  larger  amount  of 
cooling  than  the  corresponding  two-dimensional  disturbance  terms. 
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Rccently,  Duck  (12),  considering  the  effect  of  curvature  on  the 
inviscid  axisymmetric  stability  of  a  compressible  boundary  layer, 
determined  that  it  has  a  profoundly  stabilizing  effect.  Considering 
acisymmctric,  inviscid  disturbances  in  the  supersonic  boundary  layer  formed 
on  a  thin  straight  circular  cylinder  with  adiabatic  wall  conditions,  he 
determined  that  the  effect  of  curvature  is  to  completely  stabilize  the 
''first  mode''  of  inviscid  instability  at  a  comparatively  short  distance 
down  the  axis  of  the  cylinder  and  causes  the  ''second  mode''  of  instability 
to  undergo  substantial  maximum  growth  rate  reductions  at  increasingly 
further  downstream  locations  along  the  cylinder  axis.  All  the  numerical 
evidence  though,  is  that  for  large  distances  downstream  the  ''second  mode' 
is  still  present,  having  a  very  small,  but  measurable  growth  rate.  Duck 
also  observered  that  the  inclusion  of  curvature  terms  causes  the 
generalised  inflexion  condition  described  above  to  be  altered  and  he 
derived  a  modified  (or  ’’doubly  generalised'')  inflexion  condition, 
involving  the  radius  of  curvature. 

The  work  of  Duck  (12)  was  extended  to  include  non-ax i symme t r i c 
disturbance  terms  by  Duck  and  Shaw  (13)  (here-after  refered  to  as  DS) 
and  they  applied  their  techniques  to  a  different  axi symmet ric  configuration, 
namely  that  of  a  sharp  cone.  DS  generalised  the  doubly  generalised 
inflexion  condition  further,  to  give  a  condition  for  the  existence  of  so 
called  inviscid,  neutral  non-axi symmet ric  modes  and  termed  it  the 
'triply  generalised  inflexion  condition’’. 

In  this  paper  we  consider  the  combined  effects  that  curvature 
and  wall  cooling  (or  heating)  have  on  the  stability  of  the  compressible 
boundary  layer  associated  with  an  axisymmetric  body  (principly  a  thin, 
straight  circular  cylinder).  We  are  interested  in  both  axisymmetric  and 
non- ax i symme t r ic ,  inviscid  disturbances.  The  theory  for  the  form  of  the 
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boundary  layer  and  the  resultant  disturbance  equations  has  already  been 
presented  by  Duck  (12)  and  in  DS ,  to  which  we  refer  the  reader  to  for 
details,  but  we  shall  begin  this  paper  with  a  brief  account  of  the  relevant 
equations  to  the  present  problem.  It  should  be  noted  that  although  the  work 
in  DS  was  for  a  sharp  cone,  the  necessary  equations  for  a  circular  cylinder 
for  this  problem  are  easily  obtained  by  setting  the  slope  parameter,  X 
in  DS  to  zero. 


§2  Formulation  of  Problem 

The  general  layout  of  the  problem  is  shown  in  Figure  1.  The  z* 

axis  lies  along  the  cylinder  axis,  r*  denotes  the  radial  coordinate. 

and  8  the  azimuthal  coordinate.  We  have  a  supersonic  flow  incident 

normal  on  the  end-face  of  the  cylinder  (which  has  radius  a*),  where 

* 

H*  represents  the  free-stream  Mach  number,  U„  the  free-stream  velocity 

★  *  * 

(in  the  axial  direction)  and  p„,  p„  and  TM  represents  the  free-stream 
density,  first  coefficient  of  viscosity  and  temperature  respectively.  We 

define  the  Reynolds  number  to  be 

*  *  * 

Re  =  ik-{| — ,  (2.1) 

and  this  will  be  assumed  to  be  large  throughout.  Subscript  °°  denotes 
free-stream  conditions. 

For  the  resultant  boundary  layer  formed  on  the  cylinder  we  make  the 
following  assumpt ions : - 

Firstly  following  previous  authors  we  shall  assume  an  absence 
of  any  shock  wave  in  the  basic  flow.  We  feel  justified  in  making  this 
assumption,  as  the  chosen  'thinness'  of  the  boundary  layer  is  such  that  the 
shock  wave  will  be  sufficiently  far  away  from  the  boundary  layer  to  have 
negligible  effect.  Secondly  we  shall  assume  the  linear  Chapman  viscosity 
law,  where 

4*  =  C  T*.  (2.2) 

and  C  is  assumed  to  be  constant.  Since  the  flow  is  one  of  high 
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Reynolds  number,  wc  make  the  steady  bounda ry - 1  aye r  approximation.  A 
fundamental  component  of  this  work  is  the  inclusion  of  curvature  terms 
into  the  governing  equations.  To  achieve  this  we  assume  that  the  body 
radius  is  generally  of  the  same  order  as  the  boundary- l ayer  thickness  (a 
similar  approach  was  adopted  by  Seban  and  Bond  (14),  Duck  and  Bodonyi 
(15),  Duck  (12)  and  in  DS),  except  at  the  tip  of  the  cylinder  where  we 
shall  assume  the  boundary  layer  has  zero  thickness  (as  assumed  by  Seban  and 
Bond  (14)).  The  final  assumption  we  shall  make  is  that  of  a  perfect  gas. 

Making  use  of  the  above  mentioned  assumptions,  in  DS.  we  obtained  the 
following  system  of  non-dimensional,  leading  order  equations  for  the  basic 


(2.3) 


(2.4) 


^3  +  v3  ^3  =  I  L  [  r  T 
Jr  dz  r  dr  1 


(2.5) 


2  f  dv-i 


v1-  +  V3-  =  T2(Y-l)M^(-l) 


T  d  r  r  dT  i 

+  7  JF  [  a  T  37  J. 


(2.6) 


where  U„(Re' ^vj ,V2,V3)  represent  the  velocity  components  in  the 
(r*=a*r,  0,  z*  =  Re  a*C'^z)  directions  respectively,  p*R*T*P 
(where  R  is  the  gas  constant)  and  T„T  represent  the  pressure  and 

temperature  terms  respectively,  y  denotes  the  ratio  of  specific  heats 

p  CD  * 

and  o  =  — j— K  is  the  Prandtl  number,  with  K  being  the  coefficient  of 
K 

heat  conductivity;  in  (2.3)-(2.6)  the  basic  flow  is  assumed  independent  of 
9,  and  has  no  azimuthal  velocity  component  (i.e.  V£  =  0). 

The  boundary  conditions  for  the  problem  are 


V1  =  v3  =  0  on  r=l . 

V3  -*  I ,  T  -*  1  as  r-»«>. 


(2.7) 


together  with  a  wall  temperature  condition  which  for  heated  or  cooled  walls 
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has  the  form 

T  =  Tw  on  r=l .  (2.8) 

The  system  of  equations  (2.3)-(2.6)  was  solved  using  a  straight¬ 
forward  Crank-Nicol son  scheme,  as  described  in  DS,  subject  to  the  boundary 
conditions  (2.7)  and  (2.8). 

To  study  the  stability  of  this  basic  flow,  we  now  investigate  the 
effect  of  small  amplitude  disturbances.  We  shall  assume  that  the 
disturbance  wavelength  is  generally  comparable  to  the  boundary- 1 ayer 
thickness  and  therefore  also  of  the  body  radius  ( 0 ( a  * ) ) .  in  which  case 
the  parallel  flow  approximation  is  completely  vindicated. 

At  a  fixed  z-station  we  express  the  flow  parameters  of  velocity, 
pressure,  temperature  and  density  as  the  sum  of  a  mean  flow  term  plus  a 
small,  first  order  disturbance  term,  for  example, 

v  j  =  5aU„  <p(  r)  E  +  0(52) , 

v3  =  U^twQfr)  +  5  v3( r )E 1  +  0(82). 

T*  =  T*  [T0( r )  +  5  T ( r)E]  +  0(52), 

P*  =  pVt!  [  l+5p( r)E]  +  0(52 ) ,  (2.9) 

where 

E  =  exp[  ^  (z-ct)  +  in  0),  (2.10) 

T0(r)  =  T  (r.z) 

w0( r )  =  v3(r-z)  (2.11) 

and 

c  =  zi,  T)  =  (r-l)/C,  Z  =  z*/a* .  (2.12) 

and  5  is  the  perturbation  parameter  and  as  such  is  assumed  to  be 
dimini  shingly  smal 1 . 

In  this  paper  we  focus  our  attention  on  temporal  stability  for  which 
the  growth  rate  is  acj,  where  a  represents  the  spatial  wavenumber 
and  Cj  is  the  imaginary  part  of  the  wavespeed.  If  Cj  >  0  we  then 
have  growing  disturbances,  if  Cj  =  0  the  disturbance  is  neutral  and  if 
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Cj  <  0  the  disturbance  decays. 

Substitution  of  the  flow  parameters  into  the  full  system  of  equations 
of  continuity,  momenta,  energy  and  state,  discarding  0(52)  terms  and 

all  but  the  largest  terms  in  Re,  we  obtain  a  sixth  order  system,  which 

after  some  algebra  (as  described  in  DS)  can  be  reduced  to  the  following 


system 


ia2(w0-c) 


<p  =  -  pn  . 

yM2* 


(2.13) 


whe  re 


0  =  (WQ-09T,  -  [Won  - 
YM2„  1+tiC 


*  =  To  [1  +  - -1——  ■  ]  -  m1(w0-c)2. 

a2(l+ri02 


The  boundary  conditions  are 
<p  =  p^  =  0  on  q=0, 


(2  14) 


(2.15) 


(2.16) 


where 


<P  ~  <P~  M1!) 


P  -  P»Kn(Tl) 


as  t)  -»  », 


T)  =  ±  0t[  I -m3(  I -c  )2 1*  (£  +  *1). 


(2.17) 


(2.18) 


It  should  be  noted  that  the  argument  of  the  modified  Bessel  function 
(i.e.  appropriate  sign  in  (2.18))  Kn(q)  is  chosen  to  ensure  boundedness 
in  the  farfield. 


The  system  of  equations  (2.13),  (2.14)  was  solved  using  a  Runge- 
Kutta  scheme  for  the  eigenvalue  c,  as  described  in  DS,  subject  to  the 
boundary  conditions  (2.16),  (2.17). 
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§3  Numerical  Results 

Most  of  this  section  will  be  devoted  to  the  effect  of  wall-cooling 
on  the  stability  of  compressible  boundary  layers  and  its  interaction  with 
curvature,  as  these  conditions  are  likely  to  be  of  interest  in  important 
applications.  However,  at  the  end  of  the  section  we  shall  give  some  heated 
wall  results  which  exhibit  some  additional  interesting  physical  features. 

As  in  DS  we  shall  make  one  choice  of  Prandtl  number  (0=0.72), 

Mach  number  (M0O=3.8)  and  ratio  of  specific  heats  ( y=  1.4). 

3.1.  Cooled  Wall  Results 

We  begin  by  considering  the  effect  that  wall  cooling  has  on  the 
inflexion  points.  In  DS  a  condition  was  derived  for  the  existence  of 
inviscid,  neutral,  non  -  ax i symme t r i c  modes  in  the  boundary  layer  on  an 
axisymmetric  body,  which  was  termed  the  'triply  generalised  inflexion 
condition''.  However  this  condition  does  involve  the  azimuthal  and 
streamwise  wavenumbers  a  and  n,  and  so  it  is  difficult  to  forecast, 
prior  to  any  numerical  investigation,  the  existence  of  neutral  stability 
points  of  this  kind.  However  in  the  case  of  axisymmetric  disturbances, 
this  is  no  longer  the  case,  since  the  condition  reduces  to 

d  r  wq r ( r )  1=0.  (3.1) 

17  1  -T0(r)  1 

as  determined  by  Duck  (12). 

Figure  2a  shows  the  axial  variation  of  radial  position  of  the 

generalised  inflexion  points  for  the  temperatures  shown  As  in  the 
insulated  cylinder  case  (12),  the  graphs  display  two  prominent  features: 

(i)  the  inflexion  points  occur  in  pairs  and  (ii)  there  exists  a 
critical  value  of  downstream  of  which  no  such  points  exist.  The  point 
£  =  0  corresponds  to  the  tip  of  the  cylinder  and  as  such  corresponds  to 
the  planar  case  as  studied  by  Mack  (2).  We  find  that  as  the  surface 
of  the  cylinder  is  cooled,  the  lower  inflexion  point  lifts  up  off  the 
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cylinder  surface.  For  sufficient  cooling  the  lower  point  coalesces  with 
the  upper  inflexion  point  and  further  cooling  results  in  the  complete 
disappearance  of  the  inflexion  points.  Therefore  for  a  given  ^-station 
there  exists  a  critical  wall  temperature  below  which  no  inflexion  points 
exist.  In  light  of  our  earlier  comments,  we  do  have  an  additional  condition 
for  the  existence  of  neutral  subsonic  disturbances,  namely  that 

1  ■  h  c  ‘  1  *  jl-  13  21 


This  has  direct  implications  on  the  first  mode  of  instability  as  this 
generally  requires  the  presence  of  an  inflexion  point  in  the  profile 
satisfying  (3.2). 

Figure  2b  shows  the  axial  variation  of  WQ(r)j)  for  t*1c  displayed 
wall  temperatures.  From  stability  theory,  unstable  subsonic  modes  only 
exist  if  a  generalised  inflexion  point  satisfying  (3.2)  occurs  somewhere 
in  the  boundary  layer.  Examination  of  the  curve  for  Tw  =  3.0  reveals  that 
subsonic  generalised  inflexion  points  only  occur  for  0  <  £  <  0.0795; 
consequently  the  mode  has  completely  disappeared  before  the  generalised 
inflexion  points  have  merged.  For  a  wall  temperature  of  Tw  =  2.0  the 
generalised  inflexion  points  are  all  now  supersonic  in  nature  (c  <  1  -  , 

implying  for  this  and  all  cooler  wall  temperatures  the  eradication  of  Mode  I 
instabilities. 

We  now  turn  our  attention  to  the  eigenvalue  problem  for  both 
axisymmetric  and  non- ax i symme t r i c  disturbances.  We  shall  present  only 
results  for  unstable  modes  and  all  plots  are  for  cj  (where  c  =  cr  +  icj). 
The  first  set  of  results  we  present  corresponds  to  the  tip  of  the  cylinder 
and  as  such  are  comparable  to  the  planar  results  as  obtained  by  Mack  (2,  10, 
11).  Figure  2c  displays  distributions  for  the  first  mode  of 
instability.  We  observe  that  all  the  modes  originate  as  sonic  (i.e.  with 
c  =  l  -  |Jj-)  neutral  modes  and  terminate  as  subsonic,  generalised 
inflexional  modes.  We  note  that  as  we  cool  the  surface  of  the  cylinder 
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this  mode  undergoes  stabilization,  until,  with  sufficient  cooling  it  becomes 
completely  stable,  thus  verifying  Mack’s  observations  (and  those  of  Lees 
(3),  Van  Driest  (4)  and  Dunn  and  Lin  (5)).  From  our  inflexion  point 
results,  this  is  to  be  expected,  as  we  observed  that  for  cool  enough 
wall  conditions  w0  i )  droPs  below  1  -  ^  and  we  no  longer  have  the 
necessary  conditions  for  a  subsonic  generalised  inflexional  mode. 

Figure  2d  shows  the  second  mode  of  instability  at  the  tip  of  the 
cylinder.  All  the  modes  originate  as  subsonic  generalised  inflexional 
modes  with  the  special  case  c(  =  cr  =  wg(T)j))  =  1,  which  therefore 
corresponds  to  a  generalised  inflexion  point  in  the  frecstream,  and 
terminate  as  neutral  modes.  Depending  on  whether  or  not  a  subsonic 
generalised  inflexional  mode  for  the  given  wall  conditions  exists,  this 
neutral  mode  is  either  subsonic  inflexional  or  supersonic,  in  nature, 
and  the  mode  may  continue  as  a  decaying  mode  (Cj<0).  Examination  of  our 
results  reveals  that  as  we  cool  the  cylinder  surface,  the  maximum  value 
of  Cj  increases  slightly  and  then  decreases  but  the  important  product 
acj  ,  actually  increases  with  wall  cooling.  Thus  we  deduce  that  wall 
cooling  destabilises  the  second  mode  of  instability,  in  line  with  Mack's 
(2,  10,  11)  observations. 

The  next  set  of  results  we  present  corresponds  to  a  relatively  small 
distance  from  the  tip  along  the  cylinder  in  the  axial  directir  ,  at  the 
location  £  =  0.05,  and  for  an  axisymmetric  mode  (ie,  n=0) .  Figure  2e 
displays  the  Mod  I  distributions  and  we  note  that  with  sufficient  cooling, 
it  is  again  possible  to  completely  stabilize  this  mode.  We  observe  that  the 
neutral  mode  at  which  this  instability  originates  now  occurs  for  a 
slightly  greater  than  zero  and  it  is  very  slightly  supersonic  in  nature. 

When  compared  with  the  corresponding  planar  results  we  note  that  curvature 
(as  noted  by  Duck  (12)  and  DS)  has  a  profoundly  stabilizing  effect.  Even 
though  this  station  is  only  a  relatively  short  distance  along  the  cylinder 


from  the  tip,  curvature  has  reduced  the  value  of  maximum  cj  for  the 
Tw  =  3.0  curve  by  a  factor  of  4,  while  for  the  Tw  =  2.8  curve  this 
factor  is  almost  5.  Curvature  results  in  the  mode  requiring  less  cooling 
to  completely  stabilize  it.  Figure  2f  displays  the  axisymmetric  mode  II 
instability  at  this  axial  location.  Again  we  note  that  curvature  has  had  a 
stabilizing  effect  on  the  instability,  but  cooling  causes  the  mode  to  become 
more  unstable,  in  line  with  the  planar  results  described  previously. 

The  next  set  of  results  we  present  corresponds  to  the  £=0.05  location, 
for  azimuthal  wavenumber  n=l .  Figures  2g,  2h  display  the  mode  I  and  II 
instabilities,  respectively.  We  observe  that  in  this  case  the  mode  I 
instability  is  substantially  more  unstable  than  the  axisymmetric  case.  We 
find  again,  however,  with  sufficient  cooling  we  can  completely  stabilize 
this  mode,  although  the  mode  does  persist  for  cooler  wall  conditions.  The 
mode  II  instability  has  the  same  qualitative  features  as  the  axisymmetric 
case,  although  it  is  slightly  more  stable  than  the  axisymmetric  case. 

In  DS  it  was  noted  that  near  the  cone  tip,  for  non-axi symme t r i c 
disturbance  terms,  a  third  mode  of  instability  is  seen  to  develop.  This  new 
mode  originates  as  a  neutral  mode  at  a=0,  with  Cj*0,  terminates  as  a 
supersonic  neutral  mode  and  was  termed  in  DS,  mode  IA.  As  has  already 
been  noted,  the  work  of  DS  was  for  a  cone  with  adiabatic  wall  conditions 
and  it  was  found  that  this  new  mode  had  already  amalgamated  with  mode  I 
instability  at  the  £=0.05  location  for  n=l .  We  find  however  that  for  a 
cooled  cylinder  at  this  axial  location  and  n=l ,  the  mode  IA  instability 
is  still  distinct  (Figure  2i)  over  the  range  of  Tw  shown;  a  possible 
partial  explanation  is  provided  in  Section  6.  We  observe  that  wall  cooling 
causes  mode  IA  instability  to  become  less  unstable  and  with  sufficient 
wall  cooling  we  can  completely  stabilize  it.  Comparing  with  the  mode  I 
instability  we  note  that  the  maximum  value  of  cj  for  the  mode  IA 
instability  is  an  order  of  magnitude  larger,  but  due  to  the  smallness  of  a. 


the  mode  I  growth  rates  are  generally  larger.  One  further  observation, 
though,  is  that  mode  IA  instability  persists  for  cooler  wall  temperatures 
than  does  the  mode  I  instability. 

We  now  consider  the  situation  for  an  azimuthal  wavenumber  of  n=3.  at 
this  axial  location.  We  observe  that  for  this  case  the  mode  I  and  IA 
instabilities  have  now  amalgamated.  The  new  combined  mode  originates  as  a 
neutral  mode  (but  with  c{*0  for  a=0)  and  terminates  as  a  subsonic 
generalised  inflexional  mode  (Figure  2j).  Again  we  note  that  sufficient 
wall  cooling  can  completely  stabilize  this  mode,  although  it  is  found  that 
this  mode  persists  for  cooler  wall  temperatures  as  n  increases.  Comparing 
with  the  n=l  results  we  find  that  the  increase  in  n  has  also  caused  the 
mode  to  become  slightly  less  unstable.  The  mode  II  instability  (Figure 
2k),  again,  has  the  same  qualitative  features,  but  has  also  undergone 
stabilization  with  the  increase  in  n.  However,  it  is  found  that  cooling 
has  the  more  dominant  destabilizing  effect  here. 

The  next  set  of  results  we  present  correspond  to  an  azimuthal  wavenumber 
of  n=5,  at  this  axial  location.  We  observe  the  same  qualitative  features 
as  the  n=3  results  for  both  modes  (Figures  21,  2m),  although  we  note  that 
we  have  had  a  more  significant  stabilizing  effect  due  to  the  increase  in  n. 
The  mode  I  instability  is  now  completely  stabilized  for  higher  wall 
temperatures,  while  again,  it  is  observed  that  cooling  causes  an  even  more 
marked  destabilizing  effect  for  the  mode  II  instability. 

We  find  that  as  the  azimuthal  wavenumber  n  is  further  increased,  both 
the  mode  I  and  II  instabilities  undergo  additional  stabilization, 
although  cooling  maintains  a  destabilizing  effect  on  Mode  II.  For  large 
n  it  is  found  that  Mode  I  instability  is  completely  stabilized  while  the 
mode  II  instability  persists,  but  with  much  diminished  growth  rates.  In 
the  next  section  we  shall  consider  the  asymptotic  structure  of  the 
disturbance  equations  in  this  limit. 
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The  final  set  of  results  we  shall  present  in  this  sub-section  corresponds 
to  the  axial  location  ^=0.5  and  for  an  azimuthal  wavenumber  of  n=1 .  We 
find  that  at  this  distance  along  the  cylinder  the  mode  1  and  1A 
instabilities  have  now  combined  (Figure  2n).  Comparison  with  the  n=l 
results  at  the  £=0.05  location,  we  note  that  growth  rates  have  been 
reduced  due  to  the  stabilizing  effect  of  curvature.  It  is  found,  however, 
that  the  combined  mode  prevails  for  cooler  wall  conditions,  but  again  is 
completely  stabilized  with  sufficient  wall  cooling.  Figure  2o  displays  the 
mode  II  instability.  Curvature  has  resulted  in  the  growth  rates  being 
reduced,  but  we  observe  that  cooling  has  a  more  profound  destabilizing 
effect  here. 

As  we  move  further  along  the  cylinder  in  the  axial  direction,  the  mode 
I  and  II  instabilities  undergo  further  growth-rate  reductions.  For  a 
given  wall  temperature,  Tw,  azimuthal  wavenumber  n,  axial  wavenumber  a, 
there  exists  a  critical  value  of  £,  beyond  which  no  triply  generalised 
inflexion  points  (see  DS  for  details)  occur.  Consequently,  we  expect  the 
mode  I  instability  to  have  disappeared  for  axial  distances  larger  than 
this  critical  £,  which  is  bourne  out  by  our  numerical  results.  It  is 
found  however  that  due  to  the  special  nature  of  the  inflexional  point  within 
the  freestream,  the  mode  II  instability  still  originates  with  c=l ,  even 
for  axial  distances  greater  than  this  critical  £,  although  the  mode  has 
much  reduced  growth  rates.  In  the  light  of  this,  in  Section  5  we  shall 
consider  the  form  of  the  disturbance  equations  in  the  limit  of  large  £. 

3.2  Heated  Vail  Results 

We  begin  by  considering  the  effect  that  wall  heating  has  on  the 
inflexion  points..  We  restrict  our  study  to  the  case  of  axisymmetric 
disturbances,  and  consider  the  effect  wall  heating  has  on  condition  (3.1). 
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Figure  3a  shows  the  axial  variation  of  (radial)  position  of  the 
generalised  inflexion  points  for  the  temperatures  shown.  We  observe  again 
the  features  seen  in  Figure  2a.  Close  to  the  cylinder  tip.  however,  we  find 
that  for  a  small  axial  distance  measured  from  the  tip,  we  no  longer  have  any 
lower  inflexional  points.  As  we  heat  the  surface  of  the  cylinder  further, 
this  axial  distance  is  found  to  increase.  It  is  also  observed  that  wall 
heating  causes  the  critical  value  of  £,  beyond  which  no  inflexional  points 
exist,  to  increase.  For  Tw=4.5,  the  critical  value  of  £  is  about 
0.216,  while  for  Tw=6.0,  we  have  had  a  substantial  increase  to  a  value  of 
C  =  0.423.  This  will  have  direct  implications  on  the  first  mode  of 
instability  which  we  expect  to  persist  for  larger  axial  distances  as  the 
wall  is  heated.  These  effects  are  in  many  ways  to  be  expected,  being  the 
converse  of  the  cooling  observations  described  earlier. 

Figure  3b  shows  the  axial  variation  of  wQ(pj)  for  the  displayed  wall 
temperatures.  The  most  marked  feature  of  these  curves  is  that  as  we  heat 
the  cylinder  surface,  the  lower  inflexional  point  becomes  subsonic  past  a 
critical  value  of  £,  which  is  temperature  dependent,  i.e.  for  axial 
distances  greater  than  this  critical  value  of  £  but  upstream  of  the  station 
beyond  which  no  inflexional  points  occur,  both  generalised  inflexional 
points  are  now  subsonic  in  nature.  The  first  wall  temperature  for  which  we 
observe  a  lower,  subsonic  inflexional  point  is  for  a  wall  temperature  of 
about  Tw=4.5.  It  is  found,  however,  that  the  critical  value  of  C  here,  is 
very  close  to  the  stations  where  the  generalised  inflexional  points  coalesce. 
For  Tw=5.0  we  observe  that  for  the  range  0.2635  <  C  <  0.2720  we  have  two 
subsonic  inflexional  modes,  while  for  the  hotter  wall  temperature  of  Tw=6.0 
we  have  the  larger  range  0.363  <  C  <  0.423  for  which  both  generalised 
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inflcxional  points  arc  subsonic.  In  these  ^-ranges  we  expect  the  appearance 
of  two  subsonic  generalised  inflexional  modes  to  have  a  significant  effect 
on  the  physics  of  the  problem. 

We  now  present  eigenvalue  results  for  axisymmetric  disturbances  only. 

As  in  the  case  of  cooled  wall  conditions  we  focus  attention  on  unstable 
modes.  We  begin  by  considering  the  effect  wall  heating  has  on  the 
mode  I  and  II  instabilities  for  a  ^-station  close  to  the  cylinder  tip 
(£=0.05)  and  consequently  the  lower  inflexional  point  is  still  supersonic 
in  nature.  Figure  3c  shows  the  mode  I  instability  for  the  temperatures 
shown.  We  observe  that  all  the  modes  originate  as  neutral  modes  at  a  value 
of  a  slightly  greater  than  zero,  and  are  very  slightly  supersonic  in 
nature.  We  find  as  we  heat  the  wall  this  neutral  mode  tends  towards  the 
sonic  value.  All  the  modes  terminate  as  subsonic  generalised  inflexional 
modes,  continuing  on  as  stable  modes  ( c j <0 )  for  larger  values  of  a. 

These  observations  are  similar  to  the  results  obtained  for  the  axisymmetric 
cooled  wall  case  for  this  ^-station.  We  find,  as  we  would  expect,  heating 
the  surface  of  the  cylinder  causes  the  mode  I  instability  to  become  more 
unstable  -  converse  to  the  effect  of  cooling  on  this  mode. 

Figure  3d  displays  the  mode  II  instabilities  for  the  £=0.05  station, 
at  the  temperatures  shown.  We  find  all  the  modes  originate  as  subsonic 
generalised  inflexional  modes,  rise  to  a  maximum  and  terminate  as  subsonic 
generalised  inflexional  modes  (which  then  continue  in  all  the  cases 
presented  as  stable  modes).  Heating  the  cylinder  wall  causes  the  mode  II 
instability  to  become  less  unstable  and  the  numerical  evidence  suggests  that 
with  sufficient  heating  we  can  completely  stabilize  this  mode. 

We  new  consider  the  effect  that  the  lower  inflexional  point  becoming 
subsonic  has  on  the  mode  I  and  II  instabilities.  Figure  3e  displays 
the  mode  I  instability  for  a  wall  temperature  of  Tw=5.0  and  for  the 
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£-stations  as  marked.  For  the  £=0.26  station  (where  the  lower  inflexional 
point  is  still  supersonic)  the  mode  originates  as  a  very  slightly  supersonic 
mode  and  terminates  as  a  subsonic  generalised  inflexional  mode.  For  the 
£=0.264  station  the  lower  inflexional  point  has  now  become  subsonic  in 
nature,  and  we  find  the  mode  I  instability  now  originates  as  this  lower 
subsonic  generalised  inflexional  mode.  Consequently  the  value  of  a  for 
the  neutral  mode  has  increased  correspondingly.  As  before,  the  mode  I 
instability  terminates  as  the  upper  generalised  inflexional  mode  which  is  of 
course  subsonic,  as  well.  From  the  inflexional  point  curves  we  know  that  as 
we  move  upstream  the  inflexional  points  move  closer  together,  eventually 
coalescing  and  this  is  reflected  in  the  new  form  of  the  mode  I 
instabilities.  For  the  £=0.27  station  the  mode  I  instability  occurs 
over  a  much  smaller  a-range  and  the  growth  rates  are  greatly  reduced. 

Figure  3f  displays  the  mode  I  instability  for  a  wall  temperature  of 
Tw=6.0  and  the  displayed  £-stations.  Again  we  note  that  as  the  lower 
inflexional  point  becomes  subsonic  the  neutral  point  at  which  the 
instability  originates  transforms  from  being  very  slightly  supersonic  in 
nature,  to  this  inflexional  mode.  Movement  upstream  causes  the  a-ranges 
and  growth  rates  to  be  diminished,  but  the  reduction  is  less  marked  (in 
comparison  with  the  Tw=5.0  results)  due  to  the  destabilizing  effect 
brought  on  by  wall  heating. 

The  appearance  of  a  second  subsonic  generalised  inflexional  mode  is 
found  to  have  no  effect  on  the  mode  II  instability  as  it  always  terminates 
as  the  upper  generalised  inflexional  mode.  It  is  found,  however,  that  a 
third  mode  of  instability  exists,  originating  as  the  lower  generalised 
inflexional  mode  and  terminating  as  a  slightly  supersonic  neutral  mode. 

This  new  mode,  which  we  shall  term  Mode  1 1 A .  occurs  for  values  of  a 
greater  than  the  value  of  o  for  which  the  mode  II  instability  terminates. 

It  appears  that  the  mode  II  instability  continues  as  a  stable  mode  and 
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thcn  becomes  unstable  again  at  the  lower  generalised  inflexional  point. 

The  growth  rates  of  the  mode  1 1 A  arc  found  to  be  very  small.  For  a  wall 
temperature  of  Tw=5.0  and  the  station  £=0.27  the  growth  rates  are  in  the 
order  lO'll ,  while  for  Tw=6.0  and  £=0.364  the  growth  rates  are  of  the 
orde r  of  10"  *  3  _  io*  14 

We  now  move  on  to  consider  the  form  of  the  disturbance  equations  in 
the  limit  of  large  azimuthal  wavenumbers. 


§4.  Disturbance  Equations  for  Large  n 

In  this  section  we  consider  the  form  of  the  disturbance  equations  in  the 
asymptotic  limit  of  large  n,  guided  by  our  numerical  results. 

4.1  Formulation  of  the  Problem 

Consider  the  pressure  disturbance  equation  as  derived  in  DS 


whe  re 


^  !_  r  IoJji  1  +  r  _  W0T1]  =  OP. 

a/-  dq  L  wq-c  J  L  l+q£  /T'J  a^fwg-c) 

*  -  T°  C 1  +  527T&  3  -  M- 


(4.1) 


(4.2) 


In  the  limit  of  large  n  our  numerics  suggest  that  the  corresponding 
streamwise  wavenumbers  for  the  instability  also  increase,  and  that 
<t>(q=0)  — »  0.  Equation  (4.2)  suggests  we  must  have 


a  =  an,  a  =  0(1) . 


(4.3) 


A  (sensible)  balancing  of  terms  in  equation  (4.1)  suggests  asymptotic 
expansions  of  the  form 

c  =  eg  +  n*2^3  cj  +  0(n'4^3)  +  ...  , 


4>  =  <t>o  +  n*2/3  Oj  +  0(n*4/3)  +  ...  . 


(4.4) 
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The  leading  order  pressure  equation  has  the  form 


iVl  -  ¥5  *0  p  -  o. 


(4.5) 


whe  re 


*0  '  To  [  '  *  1  -  »l  <»o-c0)2- 


(4.6) 


This  has  a  solution  of  the  (bounded)  WKB  type  (assuming 

Im{(-*0>*}  >  0 


p  =  -J|}r-  {  ‘  T§r  J  [  -  *o  ]  *  dr\  } 


(4.7) 


We  can  re-write  the  expansion  for  <D  in  the  neighbourhood  of  q  =  0  as 


d>  =  d>0(0)  +  n‘2/3  [  <J>0  (0)  q  +  4>,  (0)  ]  +  0(n"4/3) , 


(4.8) 


whe  re 


«0(0)  -  T.  [  I  +  -g  J  -  M2  c2. 


(4.9) 


0,(0)  =  -  2  c0  c,  M‘, 


(4.10) 


and  for  insulated  wall  conditions 


v°> 


2  C3  Tw  t  u2  ' 
_ - ,-W  +  2  M  w_ 


0(0)  c0. 


(4.11) 


while  for  heated  or  cooled  wall  conditions  we  have 


V°>  -  W0)  t 1  *  &  1  -  Hr1*  * 2  H-  *o<0)  co’ 


(4.12) 


where  Tw  represents  the  wall  temperature. 
If  $g(0)  =  0,  then  we  must  have 


T  i  [  1  +  £  ]* 

w  L  TJ  -I 


(4.13) 


(which  is  clearly  real).  We  now  seek  to  determine  the  first  order 
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correction  term  to  c,  namely  C[.  For  this  we  must  study  the  wall  layer, 
which  by  a  balancing  of  terms  of  (4.2)  suggests  is  given  by 


q  =  q  n"2/3 ,  n  =  0(1). 


(4.14) 


The  (bounded)  solution  for  P  in  the  wall  layer  is  given  by 


P  =  P0  Ai  (q) , 


(4.15) 


where  Pq  is  a  constant  and  q  is  given  by 


-  f  S2  *1/3  r  '  - 

* =  77!-  [V°^V°>] 

l  Tw[<Dn(0)J  J 


(4.16) 


The  boundary  condition  at  the  wall  is 


P*1  |q=0  =  °* 


(4.17) 


which  we  transform  to  q-space,  giving  the  equation 


P-  [  I1'3  q  -  B  ]  .  0. 


(4.18) 


whe  re 


-2 

l  =  - “ - T  .  B  =  $1  (0) 

Tw[vo)r 


(4.19) 


Now  since  the  solution  in  this  region  to  the  pressure  is  given  by  Airy's 
function  (4.15)  we  have  that 


Ai  fj  [tJ  =  bX  1/3]  =  0. 


(4.20) 


Transforming  back  to  q-space  we  have  the  result 

52  1,/3 

Yn  =  -  ?  -2  [  "  *1<0>  ]• 

Tw[V0) 


(4.21) 


where  yn  (where  n  =  1,2,...)  represent  the  solutions  of  the  equation 


Ai '(-  Yn)  =  0. 


(4.22) 
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Substituting  either  (4.11)  or  (4.12)  ;i nil  equation  (4.10)  into  equation 
(4.21)  yields  the  first  order  correction  term  lor  C)  namely 


j  1/3 

1  w 


C1 


{  2  K 


V'» 


C() 


1  ) 


2/3 


2  a 


TTT 


M"  c. 


Yn- 


(4.23) 


for  insulated  walls,  whilst  lor  heatcd/coolcd  wall  conditions  we  have 

T»'/3  {  2  [«;  .„(»)  c„  -  ^  ]  *  T„.0,  [  1  *  £  ]  }2/3 

c,  =  - — - “ -  Yn-  (4-24> 

2  ( t2/3  M“  cn 

where  Cj  is  real  for  both  cases. 

Since  yn  is  a  solution  ol  equation  (4.22),  where  Yn  >  0.  then  we 
have  an  infinite  number  of  discrete,  real  possible  values  for  Yn  as  t*ie 
derivative  of  the  Airy  function  has  an  infinite  number  of  discrete  roots, 
confined  to  the  negative  real  axis  This  suggests  that  there  arc  an 
infinite  number  of  modes. 

We  shall  now  compare  these  results  «nli  numerically  determined  results 
for  large  values  of  n. 


4.2  Numerical  Results 

All  the  results  presented  in  this  section  arc  for  a  frccstrcam  Mach 
number  of  3.8  and  at  the  point  C,  =  0.2  along  the  cylinder.  We  only 
present  results  for  adiabatic  wall  conditions. 

Firstly  let  us  consider  the  asymptotic  expansion  for  c  in  the  limit 
of  large  n.  The  leading  order  term  in  the  c  expansion  eg.  is  given  by 
equation  (4.13)  and  using  the  numerically  determined  values 


Tw  =  3.379. 
a  s  0. 1325 , 


(4.25) 


we  have 


cq  =  0.7978. 


(4.26) 
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The  first  order  correction  term,  which  for  adiabatic  wall  conditions 
is  given  by  equation  (4.23),  is  computed  to  have  the  value 

Cj  =  0.3698  yn.  (4.27) 

where  we  have  used  tt'0r|(01  =  0.1904  and  yn  are  the  zeros  of  the 
derivative  of  the  Airy  function.  The  first  six  possible  values  of  Yn  arc 
determined  from  tables  (Abramowitz  &  Stegun  (16))  and  the  corresponding 
values  of  cj  arc  shown  in  table  I 


Cl 

0.3698 

(1) 

1.2092 

(11) 

1.7750 

(Ul) 

2.5146 

(IV) 

2.9214 

(V) 

3.3282 

(VI) 

TABLE  I:  VALUES  OF  q 
-2/3 

Figure  4a  shows  a  plot  of  c(  =  eg  +  n  Cj).  as  determined 
asymptotically,  against  n  for  the  different  values  of  the  first  order 
correction  term  C|,  where  the  numbering  refers  to  the  numbering  of  the 
correction  terms  in  Table  I.  (It  should  be  noted  that  only  n  integer  has 
physical  significance,  although  Figs.  4a  and  4b  show  c  as  a  continuous 
function  of  n.)  From  hereon  we  shall  refer  to  the  different  values  of 
c  as  order  1  to  VI  inclusively  corresponding  to  the  numbering  convention  of 
the  correction  terms  in  Table  I. 

Now  as  observed  above,  our  asymptotic  analysis  suggests  the  existence 
of  an  infinite,  discrete  number  of  possible  values  for  c.  When  we  searched 
for  the  eigenvalues  numerically,  for  large  n,  we  determined  that  there 
were  indeed  many  modes.  Figure  4b  displays  two  plots  of  cr  against  n 
for  order  I  and  order  V  correction  terms.  Graph  (1)  in  each  case  represents 
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the  asymptotic  curve  and  graph  (2)  is  the  numerically  determined  curve  It 
should  be  noted  that  in  this  range  of  n  and  u,  icji  «  1 
(Cj  -  10-10  _  io-12)_  comparable  to  the  machine  accuracy  of  our 
computations.  From  our  two  sets  of  plots  we  note  that  we  have  good 
agreement  between  the  numerical  solutions  and  asymptotic  theory  for  large 
n . 

Turning  our  attention  now  to  the  form  of  the  pressure  disturbance  terms, 
as  obtained  numerically,  we  find  that  they  do  indeed  follow  the  pattern 
predicated  by  our  asymptotic  theory,  being  initially  oscillatory  in  the 
Airy  solution  region  and  then  decaying  to  zero  in  the  far  field.  It  is 
also  observed  that  increasing  the  order  of  the  correction  term  has  the 
effect  of  increasing  the  number  of  zeros  of  the  e i genso 1 u t i on .  Figure  4c 
displays  the  distributions  of  Real  { P }  for  n  =  40  corresponding  to 
the  orders  as  shown. 

Examining  the  c  expansion  again,  we  have  determined  that  both  cq 

and  cj  are  real  and  therefore  the  leading  order  imaginary  term,  Cj, 

-4/3 

is  at  most  0(n  ).  This  means  that  the  leading  order  term  in  the  growth 

rate  (a  Cj)  is  0(n  at  most.  Therefore  actual  growth  rates  will 

decrease  as  n  — ►  »  which  is  confirmed  to  be  true  by  our  numerical 
obse  rva  t i ons . 

We  now  turn  our  attention  to  the  form  of  the  disturbance  equations  in 
the  limit  £  — *  ■». 


§5  Disturbance  Equations  for  Large  £ 

In  this  section  we  consider  the  form  of  the  disturbance  equations  in 
the  far  downstream  region,  guided  by  our  numerical  observations. 

5.1  Formulation  of  the  Problem 

Consider  the  pressure  disturbance  equation,  as  presented  in  the  previous 


section,  but  written  in  terms  of  'r'  rather  than  'n’ . 
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whe  re 

*  =  To  [  1  +  ]  -  M!  <wo-<>2. 

and 


(5.1) 


(5.2) 


r  =  1  +  ri  C  •  (5.3) 

In  the  limit  (,—**>,  we  assume  a  scale  of  a  of  the  form 

a  =  a*  (5.4) 

_ ^ 

where  a  is  to  be  determined. 

Guided  by  Duck's  (12)  work  for  the  form  of  the  basic  flow  in  the 
farfield  we  define  the  (small)  parameter 

<  =  U  log  Z)‘l  =  (log  O'1  •  (5.5) 

In  the  slow  moving  viscous  region  close  to  the  wall  (namely  the  r=0(l) 
lengthscale)  we  expect  asymptotic  expansions  of  the  form 

c  =  eg  +  e  c i  +  0(e2) , 

<J>  =  $>0  +  e  4>j  +  0(e2), 

T0  =  Tw  +  <  Tj  +  0(«2) , 

w0  =  <  wq  +  0(e2) .  (5.6) 

Now,  our  numerical  observations  strongly  suggest  that  <J»q  — *  0  as  £ 
becomes  large  in  the  r  =  0(1)  region  (which  may  be  confirmed  a  posteriori) 
which  means  that  $  =  0(e).  Examination  of  the  <t>  expression  reveals  that 
this  is  only  possible  in  general  if 


o'  =  0(e-i).  (5.7) 

Therefore  we  can  redefine  the  scale  (5.4)  to  have  the  form 

a  =  a  C  (log  (5.8) 
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whc  rc 


a  =  0(1). 

To  leading  order,  equation  (5.2)  reduces  to 

♦o'T.-M^c*.  (5.9) 

but  in  view  of  our  comments  above  that  4>q  — *  0,  as  £  — *  «,  for  r  =  0(1), 
then  we  have 

=  <5  ,0> 

which  means  cq  is  real. 

At  the  first  order  in  «  equation  (5.2)  has  the  form 

<t>l  =  Tj  +  -n-j-  (2  Cq  c  j  —  2  Cq  wq)  ,  (5.11) 

while  the  0(«)  correction  to  the  pressure  equation  (5.2)  is  given  by 

siPrr  Pr  *  t 2  '°  M-  '  *<»  -T1  p  ■  0  <5  ,2) 

We  now  transform  equation  (5.12)  using  the  same  transformation  as 
applied  by  Duck  (12)  for  the  basic  flow.  Firstly  we  employ  the  transform 

7  =  lnr ,  (5.13) 


giving 

A  A  2 

e-2f  Prr  +  [  2  CQ  (c J  -  w0)  -  Tj  -  e'2T  ]  P  =  0.  (5.14) 

Secondly  we  use  the  transform 


(5.15) 


but  since  Tw  is  constant  with  respect  to  f  we  can  simplify  (5.15)  to 
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r  =  K  Tw.  (5. 16) 

where  we  take  the  constant  of  integration  to  be  zero.  Equation  (5.14)  can 
be  re-written 

PRR  +  {  &2  T*  C2  ^0  M!  (ci  -  *)  -  Tj]  e2*  T*  -  n2  T2  }  P  =  0.  (5.17) 

where  we  have  made  use  of  the  result  obtained  by  Duck  (12)  for  the  basic 
flow 

WQ  =  E.  (5.18) 

Equation  (5.17)  is  solved  numerically  to  obtain  a  value  for  cj  subject  to 
the  condition  at  the  wall 

Pff|  =0.  (5.19) 

Rk=o 

and  that  P  is  bounded  in  the  far-field.  The  second  condition  is  obtained 
by  taking  E  — limit  of  (5.17).  i.c. 

-  ji2  if  e2*  Tw  F  =  o.  (5.20) 

where  p.2  =  2  cq  M2  a2  Tw.  (5.21) 

We  find,  to  leading  order,  this  equation  has  a  decaying  solution  of  the 

form 

P  -  ii**  K  '1/4  exp  [  -  JJ-  eR  -  \  IT  Tw  ]  .  (5.22) 

5.2  Numerical  Results 

All  the  results  presented  in  this  section  correspond  to  a  free-stream 
Mach  number  of  3.8,  and  azimuthal  wavenumber  n=l . 
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From  our  numerical  observations  we  obtain  a  value  for  the  leading  order 
term  in  the  c  expansion  (5.6)  of  the  form 

cq  =  0.4617922.  (5.23) 

Using  a  fourth  order  Runge-Kutta  scheme  equation,  (5.17)  was  solved 
subject  to  conditions  (5.19)  and  (5.22)  to  determine  the  eigenvalues  Cj.  We 
find  that  for  a  given  value  of  ft  there  appears  to  be  a  large  number  of 
discrete,  real  values  for  cj.  Figure  5a  displays  a  plot  of  C|  against 
ft  corresponding  to  the  first  five  modes,  as  shown. 

As  observed  above,  our  asymptotic  analysis  implies  the  existence  of 
a  large  number  of  discrete  possible  values  for  cj,  which  in  turn  implies 
the  existence  of  a  large  number  of  discrete  values  for  cr.  When  we 
searched  for  the  eigenvalues  by  solving  the  full  system  of  equations 
numerically,  for  large  we  determined  that  there  were  indeed  many  modes 

and  we  managed  to  identify  the  first  five  modes.  Figure  5b  displays  a 
comparison  between  the  asymptotically  determined  value  of  cr  and  the 
numerically  determined  value  of  cr  against  £  corresponding  to  the  first 
mode.  We  have  relatively  good  agreement,  since  the  error  term  in  the 
asymptotic  theory  is  0(<2),  which  is  quite  large.  Therefore  the  numerical 
results  seem  to  confirm  our  asymptotic  theory. 

The  asymptotic  theory  presented  above  tells  us  nothing  about  cj  and 
consequently  reveals  no  information  about  the  growth  rate  acj;  such  an 
investigation  would  require  a  prohibitive  amount  of  algebra. 

We  now  move  on  to  consider  the  form  of  the  disturbance  equations  in  the 
limit  of  small  (large)  £  for  cooled/heated  wall  conditions  on  a  cone. 


-28- 


§6  Disturbance  Equations  for  small  £. 

In  this  section  we  consider  the  effect  wall  heating/cooling  has  on  the 
asymptotic  theory  developed  by  DS  for  the  limit  £  — *  0  (and  also  with 
some  simple  modification  the  limit  £  — ►  «  in  the  case  of  a  cone)  for  the 
description  of  mode  IA.  It  should  be  noted  that  the  analysis  in  this 
section  is  also  applicable  for  a  cone,  the  circular  cylinder  case  being 
retrieved  when  X  =  0.  We  find  that  when  we  impose  heated/cooled  wall 
conditions,  instead  of  adiabatic  conditions  as  treated  by  DS,  the  asymptotic 
theory  is  rather  different.  Here  we  consider  the  case  of  a  — »  0,  £  — *  0, 
In  DS  it  was  shown  that  the  disturbance  equation  for  <p  in  the  limit 
a  — »  0  (assuming  n  *  0)  reduces  to 

3tI  {  11  +  ^  +  ^)lt(1  +  K2  +  ^  +  Cfl  ■  WQn[1-;f  Tp2"--"12  } 

=  n2Tj»2-*  (»0-c).  (6.1) 

In  the  limit  £  — »  0,  we  assume  expansions  for  the  =  0(1)  scale  of 
the  form 

c  =  +  C2C2  +  C3C3  +  ...  ,  (6.2) 

9  =  90^)  +  C<Pl (n)  +  C292<’1)  +  C393(n)  +  ...  .  (6.3) 

w0  =  W00<T1)  +  ^Ol^)  +  C2w02(11)  +  C3w03(tO  +  •••  •  (6-4) 

TO  =  Too(H)  +  C  T01(tD  +  C2T02Ol)  +  C3T03Ol)  +  ...  .  (6.5) 

where  wgg(q)  and  Tgg(q)  represent  the  planar  values  of  the  velocity  and 
temperature  profiles  respectively,  and  wg|(q)  and  Tgi (U )  etc.  correspond 
to  the  perturbations  to  the  basic  flow  caused  through  curvature.  A  more 
complete  expansion  for  <p  is  given  in  DS. 

To  leading  order,  we  have  for  ti  =  0(1)  that 


90  =  A0  W00(T1)- 


(6.6) 
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where  Ag  represents  an  arbitrary  amplitude  parameter. 

At  the  next  order  in  £  we  obtain  the  equation  for  <pj 

w00  fltl  -  A0  w00n  C]  +  A0  Wgg2  _  wggn  <p, 

(6.7) 

+  Ag  wg,  wggn  -  Ag  wg,n  wgg  =  kj  Tgg  . 

where  we  have  used  equation  (6.6)  and  kj  is  a  constant.  Setting  q  =  0 
in  equation  (6.7)  and  assuming  1 1 =  0  we  ^ave 

-  Ag  ci  wgg^  (q=0)  =  kj  Tgg  (q=0) .  (6.8) 

The  boundary  conditions  as  q  — *  »,  must  be  compatible  with  the  form 
of  system  (6.1)  in  the  far-field  (as  presented  in  DS),  together  with  (6.6). 
Defining 

7  =  1  +  q  C  =  0(1).  (6.9) 

then  we  must  have  an  outer  solution  of  the  form 

<p°ut  =  A  f'"'1 .  (6.10) 


where 


A  =  Ag  +  (,  Ai 


+ 


and  so  we  must  also  have 

out  |  ,  ... 

^lqjq  _  -  =  ‘for  |r=1  =  -(n+1  )Ag. 


(6.11) 


(6.12) 


Substituting  this  into  (6.7)  yields 


ci 


n  Tgg  (q=0) 
w00q 


(6.13) 


In  order  to  obtain  complex  values  of  c  in  our  asymptotic  analysis  we 
must  consider  higher  orders  in  £. 

At  the  next  order  of  £  we  have  the  system 
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w00  <P2ri  +  2X  w00  90q  +  2*l  *00  9lq  +  w02  90Tj 

-  w02rj  90  +  W01  9lq  +  n2  *00  90q  "  ci  (Pit] 

+  2  wqi  n  <Pon  +  woi  <90  -  woin  9i  -  2ti  *oiq  9o 

-  2o  cj  <pon  -  ^2  <Pon  +  woo  9i  +  n  *oo  9o 

-  ci  q>o  -  woon  92  -  2 woon  9o  -  2ti  9i  woon 

-  O2  *00q  90  =  k2  T00  +  n2  T00  f  dIl 

3  o  00 

+  k 1  Tqi •  (6. 14) 

We  now  determine  the  leading  order  imaginary  component  of  the  complex 
wavespeed  c.  Since  (6.14)  contains  only  real  coefficients,  any 
imaginaries  must,  of  necessity,  only  arise  at  a  critical  point,  where 
c  =  wqq.  Since  c  =  0(£),  this  must  occur  when  t)  =  0(£).  We  therefore 
consider  a  thin  layer  relative  to  the  q  =  0(1)  scale,  namely 

y=Ti/C  =  0(l).  (6.15) 

On  this  scale,  the  expansion  for  <p  is  expected  to  develop  as 

9  -  C  $o  (y)  +  C2  (y)  +  C3  $2  (y)  +  •••  •  (6.i6) 

where  the  are  expected  to  be  normalised  in  such  a  way  as  to  be 

generally  0(1)  quantities. 

The  leading  order  equation  has  the  form 

*0  =  A0  y  *00q  (11  =  0)  (6.17) 

(where  Ag  was  introduced  in  equation  (6.6)). 

At  the  next  order  in  £  we  obtain  the  system 

(w00tj(0)  7  -  Cj]  <5iy  -  w00T1(0)  =  kj 

+  i  V2  *00tn<0)  Ag  wgg^O)  -  T00ti(0)  y.  (6.18) 

where  kj  is  a  constant  term.  This  system  cannot  be  simplified  any 
further,  as  in  DS,  since  the  condition  Wgg^^(q=0)  =0,  no  longer  holds 
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in  the  present  (i.e.  heated/cool ed  wall)  ease. 

As  in  DS,  we  assume  Aq  and  icj  to  be  real  quantities,  as  we  arc 

at  liberty  to  do  so,  and  now  consider  (where  i  denotes  an  imaginary 

component).  This  quantity  is  triggered  by  the  +irc  jump  in  the  logarithm 
across  the  critical  layer.  We  can  re-write  (6.18)  in  the  form 

[  w00q ( 0 )  7  -  cj  ]  -  $1  w00n(0)  =  R.  (6.19) 

where 

R  =  *1  +  *  y2  w00qT1(0)  A0  wQOt1(0)  -  A°  ^5% '°)tq0t1(0)  9' 

(6.20) 

and  find 

y 

$1  =  ["00ti(0)  S'  '  C(]  f  - - *  ;  •  (6.21) 

Jo  twooT)(0)  y - e 1 1 

Evaluating  this  integral,  taking  just  the  imaginaries  together  with  the 
1 imi t  as  y  — *  «  yields 

*1  -  7t  Aq  c,  [W()0n(0)  y-Cj]  |  wQQlUl(iD - TQ0t1^?^ -  \  (6.22) 

i  woo2n<°)  Too(0)  woot)(°)  ' 

=  B*  (W00t,(0)  ?  -  ci), 

where  B'  =  k  Aq  cj  {  -.QQnn^ - ^QQn ( ^ ) -  \  .  (6.23) 

’  w002t)(0)  T00(0)  woor| ( 0 )  J 

Equation  (6.22)  now  provides  a  lower  boundary  condition  for  systems  (6.7) 
and  (6.14). 

Assuming  (6.7)  contains  real  coefficients  only  (cj  is  assumed  real, 
an  assumption  that  may  be  justified  a  posteriori)  we  have 

<p[  =  B'  w00(Tl),  (6.24) 

where  B'  is  defined  in  (6.23). 
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Considcring  the  imaginary  part  of  system  (6.14)  in  the  limit  q  — *  « 


we  obtain 


+  <P,» 


(6.25) 


By  the  form  of  <p  in  the  far-ficld  (as  presented  in  DS)  we  require 


<  l„  l„ 


-  -  (n+1 )  B' . 


(6.26) 


and  therefore 


=  -  B'n. 


(6.27) 


We  now  consider  the  imaginary  part  of  (6.14)  evaluated  at  q  =  0 


■  C1  »!i,ln=0  •  c2  ’on  U  '  ,00n<0)  ^1^0  =  k2  T°o(0)  <6  28) 


Using  the  results 


^lnLo  =  B'  w00ti<°) 


=  -  B1  ci, 
2  n=0  1 


(6.29) 


we  obtain 


C-,  = 


Wn(T1=0)  TooCtH)) 


(6.30) 


Using  the  governing  equations  of  motion  we  can  simplify  (6.30)  to  the  form 


c!  =  n2  lK  TQ02(0)  wQOnn(0) 

*oon(0) 


(6.31) 


Comparing  this  result  with  the  adiabatic  result,  as  obtained  in  DS, 
we  note  that  the  first  imaginary  term  in  c  is  an  order  in  £  larger  for 
heated/cooled  wall  conditions,  implying  larger  growth  rates  in  the  present 
case.  The  ratio  of  the  leading-order  imaginary  terms  has  the  form 
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Heated/Cooled  c;  _  _ 2  WQQnr)(0) _  ^  32) 

Adiabatic  Cj  (  wOOn^^  M+2o  ( y - 1 ) n ) 

We  find,  however,  that  altering  the  wall  conditions  from  adiabatic 
conditions  as  treated  by  DS.  to  heated/cooled  wall  conditions,  has  little 
significant  effect  on  the  asymptotic  analysis  for  the  mode  1  lower  neutral 
point  for  which  a  =  0(£i). 


§7.  Conclusions 

In  this  paper  we  have  investigated  the  supersonic  boundary- 1  aye r  flow, 
and  the  inviscid  stability  thereof,  over  axisymmetric  bodies  (with 
particular  emphasis  on  a  long  thin,  straight  circular  cylinder)  subject  to 
heated  or  cooled  wall  conditions. 

The  effect  of  wall  cooling  is,  generally,  seen  to  reduce  the  importance 
of  the  ''first  mode''  of  inviscid  instability,  while  the  amplification  rates 
of  the  ''second  mode’’  of  inviscid  instability  are  increased.  Therefore  we 
have  agreement  with  the  effect  wall  cooling  has  on  planar  boundary  layers 
(Mack  (11),  for  example). 

The  converse  effect  is  observed  with  wall  heating.  The  amplification 
rates  of  the  ''first  mode''  of  instability  increases,  while  wall-heating 
causes  the  ''second  mode''  of  instability  to  be  stabilized.  Wall-heating 
may  also  cause  the  formation  of  a  second  subsonic  generalised  inflexional 
mode,  which  results  in  the  appearance  of  an  additional  mode  of  instability, 
not  found  in  adiabatic  or  cooled-wall  studies. 

Asymptotic  studies  for  large  azimuthal  wavenumbers  and  far  downstream 
distances  reveal  that  the  real  part  of  the  eigenvalue  c  is  non-unique 
in  both  limits,  suggesting  that  there  exists  an  infinite  number  of  discrete 
possible  values  for  the  real  part  of  c  in  both  these  limits,  although  the 
corresponding  values  of  the  imaginary  part  of  c  are  exceedingly  small. 
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Fig  2e  Variation  of  cj  with  a,  M„=3.8,  £=0.05,  n=0.  Mode  I 


Fig  2f  Variation  of  Cj  with  a,  M«=3.8,  £=0-05,  n=0,  Mode  II. 
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